arXiv:1501.06588vl [astro-ph.GA] 26 Jan 2015 


Draft version January 28, 2015 

Preprint typeset using style emulateapj v. 5/2/11 


THE GALEX/S^G UV-IR GOLOR-COLOR DIAGRAM: CATGHING SPIRAL GALAXIES AWAY FROM THE 

BLUE SEQUENCE 

Alexandre Y. K. Bouquin\ Armando Gil de Paz\ Samuel Boissier^'^, Juan-Carlos Munoz-Mateos^, Kartik 
Sheth^, Dennis Zaritsky®, Jarkko Laine^, Jesus Gallego^ Reynier F. Peletier^, Benjamin R. Rock^'^*^, Johan H. 

Knapen^'^*^ 

^Departamento de Astrofisica y CC. de la Atmosfera, Universidad Complutense de Madrid, E-28040 Madrid, Spain 
^Aix Marseille Universite, CNGRS, LAM (Laboratoire d’Astrophysique de Marseille) UMR 7326, 13388, Marseille, France 
^I.N.A.F., Osservatorio Astronomico di Bologna, via Ranzani 1 1-40127 Bologna, Italy 
^European Southern Observatory, Casilla 19001, Santiago 19, Chile 
^National Radio Astronomy Observatory/NAASC, 520 Edgemont Road, Charlottesville, VA 22903, USA 
^Steward Observatory, University of Arizona, 933 North Cherry Avenue, Tucson, AZ 85721, USA 
^Astronomy Division, Department of Physics, FIN-90014 University of Oulu, P.O. Box 3000, Oulu, Finland 
^Kapteyn Astronomical Institute, Postbus 800, 9700 AV Groningen, the Netherlands 
^Instituto de Astrofisica de Canarias, Via Lactea, S/N, 38205 La Laguna, Spain and 
^^Departamento de Astrofisica, Universidad de La Laguna, E-38205 La Laguna, Tenerife, Spain 

Draft version January 28, 2015 

ABSTRACT 

We obtained GALEX FUV, NUV, and Spitzer/IRAC 3.6 //m photometry for >2000 galaxies, avail- 
able for 90% of the S^G sample. We find a very tight “GALEX Blue Sequence (GBS)” in the 
{FUV— NUV) versus (AC/E—[3.6]) color-color diagram which is populated by irregular and spiral 
galaxies, and is mainly driven by changes in the formation timescale (r) and a degeneracy between 
r and dust reddening. The tightness of the GBS provides an unprecedented way of identifying star¬ 
forming galaxies and objects that are just evolving to (or from) what we call the “GALEX Green 
Valley (GGV)”. At the red end of the GBS, at {NUV— [3.6]) >5, we find a wider “GALEX Red 
Sequence (GRS)” mostly populated by E/SO galaxies that has a perpendicular slope to that of the 
GBS and of the optical red sequence. We find no such dichotomy in terms of stellar mass (measured 
by M[ 3 , 6 ]), since both massive (M* > 10^^ Mq) blue and red sequence galaxies are identified. The type 
that is proportionally more often found in the GGV are the SO-Sa’s and most of these are located in 
high-density environments. We discuss evolutionary models of galaxies that show a rapid transition 
from the blue to the red sequence on timescale of 10® years. 

Subject headings: galaxies: evolution — galaxies: photometry 


1. INTRODUCTION 

Color-magnitude and color-color diagrams are proven 
galaxy evolution diagnostic tools. Recent studies mak¬ 
ing use of GALEX and SDSS data, such as Wyder et al. 
(2007), have clearly shown the existence of a bimodal dis¬ 
tribution of galaxies in {NUV— r) vs r color-magnitude 
diagrams (CMDs), where the redder and brighter region 
is populated mainly by early-type galaxies in a “red se¬ 
quence” , and where the bluer and fainter part of the 
diagram is populated mainly by late-type galaxies in a 
“blue sequence”. 

We construct a color-magnitude and a color-color di¬ 
agram using far-ultraviolet {FUV] effective wavelength 
151.6 nm), near-ultraviolet {NUV] 226.7nm) images from 
the Galaxy Evolution Explorer {GALEX] Martin et al. 
2005), and 3.6 micron images from the images released 
by the Spitzer Survey of Stellar Structure in Galaxies 
(S^G, Sheth et al. 2010). The near-infrared (NIR) band 
is sensitive mainly to old stars, and the UV bands to 
young stars (<1 Gyr). As we show below, a combination 
of UV and NIR bands allows us to separate star-forming 
from passively-evolving galaxies, while the use of the two 
GALEX UV bands makes this analysis very sensitive to 
transitional galaxies and to the transition timescale, as 
the FUV hand is most sensitive to the presence of young 
(OB-type) stars. Our sample is based on the Spitzer Sur¬ 


vey of Stellar Structure in Galaxies (S'^G; Sheth et al. 
2010). In this study, we gather the publicly available 
GALEX data (from data release GR6/7) for these S^G 
galaxies, and measure NUV and FUV surface and asymp¬ 
totic photometry for over 2000 galaxies (^90% of the S^G 
galaxies). We then combine this with consistent surface 
and asymptotic 3.6/im photometry from S'^G. We adopt 
a value of Aq = 75kms“^ Mpc“^ for galaxies lacking a 
redshift-independent distance measurement (as in Sheth 
et al. 2010). 


2. SAMPLE 

The S'^G is a volume-limited {d < 40Mpc), magnitude- 
limited (mBcorr < 15.5 mag, corrected for inclination, 
galactic extinction, and A-correction), and size-limited 
{F> 25 >^') survey avoiding the Galactic equatorial 
planedfe] > 30°) of 2,352 galaxies, consisting of 3.6 ^m 
and 4.5 pm IRAC-band images. Galaxies covered by 
this survey account for the majority of the significantly- 
sized galaxies in our extragalactic neighborhood, includ¬ 
ing galaxies of all Hubble types and masses, in various 
environments from field to clusters of galaxies (includ¬ 
ing 187 Virgo galaxies). The HI 21 cm-redshift selec¬ 
tion, however, biases the sample against very massive 
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Fig. 1.— Left : {FUV— NUV) versus {NUV— [3.6]) color-color diagram. Morphological types are represented by color (see legend). We 
define the blue and red sequence by dividing the plot into two regions and fitting a least-squares 1-D polynomial (a line) to each region. 
The solid blue line is the weighted-by-error linear fit for the lower-left region, and the solid red line is the weighted-by-error linear fit for 
the upper-right region. The dashed lines are parallel to their solid counterpart, but translated in the {FUV— NUV) direction by ±2 sigma 
for the GBS (blue shaded region), and translated in the {NUV— [3.6]) direction by ±1 sigma for the GRS (red shaded region). The GGV 
is also shown (green shaded region). The black arrow indicates the reddening vector for an attenuation of Ay = 0.5mag. The black dots 
are measurements for the Gappellari et al. (2013) sample (section 3.2; not included in the statistics) obtained by Zaritsky et al. (2015) 
(we adopted here (/—[3.6])=0.5mag). We use a 3.6//m stellar mass-to-light ratio T 3,6 = 0.6 or log(T 3 , 6 ) =-0.22 (IMF: Chabrier; Meidt 
et al. 2014). Right: Golor-magnitude diagram using the difference between {FUV — NUV) color and the fit to the GBS in the color-color 
diagram versus Mj 3 g]. Same color-coding as in the left panel. 


ETGs (Early-Type Galaxies)^ (Lagos et al. 2014). We 
compiled a subsample of corresponding GALEX (data 
release GR6/7) tiles using GALEXView^, in both FUV 
and NUV for over 2,100 galaxies. 

FUV and NUV asymptotic magnitudes were obtained 
using the same method as in Gil de Paz et al. (2007), cor¬ 
rected for the Milky Way foreground attenuation, mea¬ 
sured for 2053 galaxies. After excluding shallow observa¬ 
tions with magnitude errors larger than unity, and bad 
photometry due to contamination from bright stars, we 
are left with 1931 galaxies. The S"^G 3.6/im photome¬ 
try has been obtained using an identical procedure but 
adapted to the specifics of IRAG (Muhoz-Mateos et al. 
2015 in prep.). All magnitudes throughout this paper are 
given in the AB system. Numerical morphological types, 
absolute R-band magnitudes and optical colors were ob¬ 
tained from HyperLeda (Paturel et al. 2003). The total 
UV coverage of our parent S^G sample is ~90%, with 
a rather uniform fraction across all properties including 
absolute R-band magnitude, distance, and morphological 
types, and therefore represents and inherits the selection 
criteria of S'^G. 

3. RESULTS AND ANALYSIS 
3.1. UV-IR Color-Color Diagram 

Using the asymptotic photometry described above, 
we construct the {FUV—NUV) versus (ARE—[3.6]) 
color-color diagram (Eig.l). We find a clear bimodal 
pattern: a very tight blue sequence, which we call the 
GALEX Blue Sequence” (GBS hereafter), and a less 
tight red sequence we call the '"GALEX Red sequence” 

^ An extension of the S^G sample, aimed to correct for the HI 
selection effect, is currently in progress. 

^ http://galex.stsci.edu/GalexView/ 


(GRS hereafter). The distribution of the GRS is per¬ 
pendicular to that seen in optical color-color diagrams 
(e.g., Strateva et al. 2001) likely due to the mass and 
metallicity dependence of the UV upturn (Boselli et al. 
2005). The more massive and metal-rich a galaxy is, the 
more important the contribution from the UV upturn 
to the UV emission and bluer the {FUV— NUV) color. 
Alternative use of this color are still being explored 
(e.g., Hernandez-Perez & Bruzual 2014; Zaritsky, Gil de 
Paz, & Bouquin 2014, 2015). 

To quantify the GBS and the GRS, we split the color- 
color plane by choosing a line that passes through co¬ 
ordinates (2,3) and (5,1) (left panel of Eig.l). We fit 
the bluer-half with a simple error-weighted linear regres¬ 
sion (the subsample size here is 1,787 galaxies), giving 
us a slope of 0.12±0.01 and a ?/-intercept at 0.16±0.01. 
The standard deviation of the sample over this range is 
rmsGBS = 0.20 while a Gaussian fit to the residuals yields 
an even smaller value of ctgbs = 0.10. 

We also fit the redder-half (144 galaxies) with a sim¬ 
ple error-weighted linear regression (in the {NUV— [3.6]) 
direction), giving us a slope of —4.21±1.62, and a 
^^sgrs = 0.45. The location of galaxies in this dia¬ 
gram is strongly correlated with the Hubble type (see 
also Eig.4). There are only two true ellipticals in the 
GBS, ESO548-023 and NGC 855. Some of the ellipti¬ 
cals found in the GBS are actually blue compact dwarfs 
(BCDs) morphologically misclassified as ellipticals. The 
latter plus the 31 later-typed ETGs (E-S0,S0) located 
in the GBS could be either in the tail of star formation 
(Bresolin 2013) or could be rejuvenated systems (Kan- 
nappan et al. 2009; Thilker et al. 2010). Finally, early- 
type spiral galaxies (SO-a, Sa, Sab) are populating the 
GBS, the GRS and the region in between (which we name 
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Fig. 2.— GBS {blue squares), GGV {green points), and GRS {red triangles) galaxies are shown, with a linear fit to the GRS galaxies 
(solid lines) and the rmsoRS = 2 (dashed lines). Left: Optical GMD {U — V) versus R-band absolute magnitude of our sample. The solid 
blue line is the red sequence fit from Stott et al. (2009) with the intercept value adjusted for comparison. Right: {B — V) versus R-band 
absolute magnitude. 


GALEX Green Valley; GGV hereafter) which is consis¬ 
tent with a scenario where these galaxies are transitional 
objects. 

3.2. UV-IR Color-Magnitude Diagram 

Fig.l (Right) shows the residuals of the UV color to the 
GBS linear fit obtained in our UV-IR color-color diagram 
(Left), {FUV- NUV) - (FUV- NUV)gbs, as a fraction 
of the absolute magnitude at 3.6/rm, Mp gj and stellar 
mass measured assuming a Tg.g of 0.6 Mq/Lq (Meidt 
et al. 2014; Rock et al. 2015 in prep.). Stellar mass is 
not the driver for the evolution of these systems in or 
out the GBS. The distribution of our samples of ETGs 
(our data and those from Zaritsky et al. 2015) is similar, 
meaning that the bias against Hl-poor ETGs in S'^G is 
not driving our results. 

3.3. Optical Color-Magnitude Diagram 

In Fig.2 we show the optical CMD, commonly used to 
study galaxy evolution (e.g., Faber et al. 2007). We show 
the position of our GBS, GRS, and GGV galaxies using 
blue, red, and green points, respectively (Fig.l). 

Despite the low number of GGV galaxies overall, they 
fall either on the optical red sequence (RS) or on the (ill- 
defined) optical green valley (GV), with only a very few 
being optically blue galaxies. This result indicate that 
GGV galaxies are evolving off the GBS much slower in 
(A[/U—[3.6]) than in {FUV— NUV), while the optical 
colors are reddened by more than [0.4 mag] 0.2 mag in 
[(U—B)] (B—V). Alternatively, the GGV galaxies are op¬ 
tically red systems that are growing a disk {towards the 
GBS), in which case the optical reddening is less than 
[0.3mag] 0.2mag in [(U—B)] (B—V). 

It is interesting to note that the majority of the galax¬ 
ies in the optical GV belong to the GBS, i.e., they are 
star-forming systems to be classified as such at UV wave¬ 
lengths. The same can be said about the blue dots found 
in Fig.2 near the optical RS. This means that either they 
are in the tail end of star formation or they are regrowing 
a disk. In particular, a disk regrowth is unlikely, given 


that the GGV galaxies are found preferentially in dense 
environments where such a process may be more diffi¬ 
cult than in the outskirts of a cluster or a lower-density 
environment (Section 4.3). 

By using the classification of interacting systems 
within S^G given by Knapen et al. (2014) we do not 
find a clear difference between the fraction of close inter¬ 
actions in the three regions of our color-color diagram. 
When splitting the sample by level of interaction and 
morphological type, low number statistics hamper deter¬ 
mining, for example, whether or not E/SO in the GGV 
are mergers or post-starbursts. 

4. DISCUSSION 
4.1. Modeling color evolution 

In Fig.3 (left), we show single stellar population (SSP) 
models (initial mass function (IMF): Ghabrier) of Char¬ 
iot and Bruzual (2007) for different metallicities and 
ages. For all cases, the (FUV — NUV) color becomes red¬ 
der until it reaches a peak (at about 1 Gyr in the case 
of Solar metallicity), and then becomes bluer again (UV- 
upturn) as the models age. For these SSP models, the 
color transition from GBS to GRS takes less than a few 
10 ® years. 

To better understand what drives the position of galax¬ 
ies in the color-color diagram, we use a grid of disk 
galaxy models from Boissier & Prantzos (2000) with 
various circular velocities V/ and spin parameters A 
(Fig.3, right panel). These models span the entire GBS 
due to their variation in star formation timescale with 
Vc and A. Besides, the reddening arrow points along 
the GBS, so internal extinction contributes little to its 
spread, since star formation timescale and dust atten¬ 
uation are nearly degenerate in these colors (Section 
4.2). The tightness of the GBS is also due to the fact 
that the effects of an episodic SFH in the FUV-NUV 
color (compared to the NUV-[3.6]) are only expected at 
very low masses (typically rather blue NUV-[3.6] col¬ 
ors). Nice examples in this regard are PGC065367 
or ES0245-007. To produce more realistic galaxy col- 
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Fig. 3.— Number density of galaxies in the {FUV— NUV) versus {NUV— [3.6]) color-color diagram, in grayscale-coded 2D bins on a 
logarithmic scale (see color bar). Each bin is 0.2 x 0.2 mag in size. Left : single stellar population (SSP) models (CB07, Chabrier IMF) for 
various metallicities are over-plotted on the number density distribution. The open squares are various age time steps for the model with 
Solar metallicity (Z=0.02), color-coded according to the following log(age) ranges: magenta 5.0 — 8.0, green 8.0 — 9.0, blue 9.0 - 10.0, and 
red 10.0 - 10.3. The age increases from left to right, from log(age) = 5.0 to 10.3, for all models. Evolved model points (log(age) > 10.0) 
are plotted in red for all models. A linear fit to these points (dash-dotted line) is traced for visual comparison with the slope of the GRS. 
The black rectangle outlines the region of the right panel. Right : The red triangles at {NUV— [3.6]) > 4 are model cluster galaxies under 
the effect of ram-pressure stripping (RPS; see Boselli et al. 2006) where the annotations indicate when the RPS started (with £ = 0 being 
the present). The model indicate that the transition timescale of galaxies from the GBS to the GRS is very short, on the order of a few 
10® years. We also show three examples, as orange, green, and red dotted lines, of BP00-I-BC03 “disk-|-bulge” models for various B/T 
ratios, where each node corresponds to a value of the B/T ratio, namely 0.0 (i.e., BPOO disk only, open circles), 0.1, 0.2, 0.3, 0.4, and 0.5, 
progressively giving bluer to redder results. The position of the maximum B/T=0.3 estimated by Laurikainen et al. (2007) for a sample 
of SO-a galaxies is shown with white dots. An offset of 0.3 magnitude was added to the {FUV— NUV) colors for all model galaxies (see 
text). Both BPOO disk models and RPS models are based on the Kroupa 93 IMF. 


ors, we add a bulge to the Boissier & Prantzos (2000) 
disks using SSP predictions with ages between 10® and 
IOIO .3 Qyj. different formation timescale r. We 
adopt a bulge model with {NUV— [3.6]) = 5.3mag, and 
{FUV— NUVjcoloi = 2.4 mag. 

Three of these “disk+bulge” models are shown in Fig.3 
(right) for various bulge-to-total mass ratios B/T. Each 
colored line represents the colors of the resulting galaxy 
model for B/T from 0.0 to 0.5 in steps of 0.1. Adding 
a BC03 SSP-based bulge to any of our disk models does 
not reproduce the reddening in (FUV — NUV) observed 
in numerous early-type spirals unless we assume a very 
high B/T (typical observed B/T of SO and SO-a galaxies 
are below 30%; Laurikainen et al. 2007). This shows 
that the position of transitional early-type spirals is not 
only due to their (red) massive bulges but also to their 
disks being red compared to the disks of galaxies in the 
GBS. This result explains the red {FUV— NUV) colors 
of some (transitional) disk galaxies. Whether those disk 
galaxies could evolve into systems with properties similar 
to those of current-day ETGs cannot be determined from 
these data alone. 

To analyze the transformation timescale of these disks 
under a more realistic scenario, we use a model that in¬ 
cludes ram-pressure stripping of the gas (RPS hereafter) 
as a potential mechanism, albeit not the only (potential) 
one, to redden a galaxy disk. This RPS model starts from 
an unperturbed disk with V). = 270kms“^ and A = 0.04 
but with additional ram-pressure to fit the properties of 
the Virgo galaxy NGG 4569 (Boselli et al. 2006). The 
points represent the epochs when the peak of RPS oc¬ 
curred, with t = 0 corresponding to the present. Both 


the disk models and the peak RPS models are shifted 
redward by 0.3 mag in {FUV— NUV) to match the ob¬ 
served color-color distribution. The need for an offset 
has been discussed in Muhoz-Mateos et al. (2011) and 
is believed to be due to uncertainties in the calibration 
of the stellar atmospheres used in the spectral synthesis. 
This latter RPS disk model fits the observed distribution 
well, especially with a moderate bulge. The transition 
timescale derived from the GBS to the GRS in this case 
would again be of the order of a few 10® years. 

4 . 2 . Internal dust reddening 

To analyze the effects of internal dust, we show the red¬ 
dening in both colors in Fig.l (left), for Ay = 0.5mag, 
assuming the attenuation law of Galzetti et al. (1994). 
This law reproduces well the relation between IR excess 
and UV color (Gil de Paz et al. 2007). Interestingly, the 
reddening arrow indicates an almost complete degener¬ 
acy in these colors between reddening and a change in 
the star formation timescale (as illustrated by the mod¬ 
els of Boissier & Prantzos 2000). We find no change 
in the distribution of offsets from the GBS as a function 
of the galaxies’ axial ratio (not shown), so the (small) 
spread of the GBS is unlikely to be due to inclination 
effects. 

4 . 3 . Environmental effects 

In Fig.4 we show, superposed on the number density 
color-color diagram, the distribution of galaxies that are 
in Virgo or in high-density regions (Laine et al. 2014). 
The relative abundance of transitional galaxies in the 
Virgo cluster is higher for the majority of morphological 
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types, suggesting that an environmental effect is present. 
There are some indications that environment plays a role 
in determining whether galaxies lie in the GBS, GGV, or 
GRS. For example, 43% (12 of 28) of Sa galaxies clas¬ 
sified as either GRS or GGV are in the Virgo cluster, 
while only 12% (16 of 137) of those classified as GBS 
are in Virgo. The relative preference of GGV and GRS 
galaxies for the Virgo environment is found for each mor¬ 
phological class separately (Table 1), although in some 
cases the results are not statistically significant. Finally, 
although the presence of early-type spirals in the GGV 
might be related to environment, this might not be sim¬ 
ply due to RPS, there might be other, more fundamental 
drivers (that also depend on environment) which deter¬ 
mine the position of galaxies in this color-color diagram 
and in the optical GMDs. One possibility might be the 
degree of uniformity of the angular momentum of ac¬ 
creted material, which would segregate galaxies by mor¬ 
phological type, mass and color in different environments 
(less uniform accretion having likely occurred in denser 
environments). 


5. CONCLUDING REMARKS 

We show that star-forming galaxies form a tight “Blue 
Sequence” in FW-to-near-infrared color-color diagrams 
(the “GALFV Blue Sequence”, or simply GBS). One 
could even define a star-forming galaxy (SFG) as an ob¬ 
ject residing in that sequence. This is a clear advantage 
over the optical or optical-IR color-magnitude and color- 
color diagrams where star-forming galaxies are mixed 
with other types of sources, such as post-starbursts or, in 
general, galaxies in the GALEX Green Valley (GGV). Sb, 
Sa, and SO galaxies are found all across the path between 
the GBS and the GRS, through the GGV. Their position 
is not only driven by their stellar mass or their bulge 
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in groups; Kawata & Mulchaey 2008), galaxy harassment 
or processes related to the degree of uniformity of accre¬ 
tion are all viable. There could be galaxies in the field 
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TABLE 1 

Counts and ratios of Virgo galaxies per morphological type bin per region 



GBS 

GGV 

GRS 

Others' 

1 

ALL 

Type'" 

Virgo^ 

all" 

ratio® 

Virgo 

all 

ratio 

Virgo 

all 

ratio 

Virgo 

all 

ratio 

Virgo 
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“"Galaxies where (FUV- NUV)< 2 cr of the CBS fit and (NUV- [3.6]) > la of the GRS fit. 

'“"Based on Hubble t-type from RC2 

'"Prom the GOLDMine database (Gavazzi et al. 2003) 

''including Virgo galaxies 
'"ratio=Virgo / all 











The GALEX/S'^G UV-IR color-color diagram 
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Fig. 4.— (FUV—NUV) versus {NUV—[S.6]) separated into nine panels by morphological type. Black circles are galaxies with no 
environmental data associated. Open circles are galaxies in the Virgo cluster. Open squares are galaxies in high-density regions as defined 
by Laine et al. (2014). The GRS, GGV, and GBS are shown in panel E. 















































































































































































